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Abstract 

The rate of accretion of matter from a solar-type star onto a primordial black hole (PBH) that passes 
through it is calculated. The probability that a PBH is captured into an orbit around a star in a 
galaxy is found. The mean lifetime of the PBH in such an orbit and the rate of orbital captures of 
PBHs in the galaxy are calculated. It is shown that this rate does not depend on the mass of the 
PBH. This mechanism cannot make an appreciable contribution to the rate of observed gamma-ray 
bursts. The density of PBHs in the galaxy can reach a critical value - the density of the mass of dark 
matter in the galaxy. 

1 INTRODUCTION 

There have been several important developments in astronomy over the past two decades associ- 
ated with the discovery of objects that are candidate black holes. Black holes with galactic masses 
(10 6 -j- 10 9 Mq), stellar masses (~ IMq), and intermediate masses (~ 1O 3 M0) have all been identified. 
Naturally, black holes with masses that are much less than a solar mass can be detected only indirectly 
(such as via gravitational microlensing) . So-called primordial black holes (PBHs) are in this class. 
They can be formed only in the early stages of evolution of the Universe, since the natural evolution 
of a star does not permit the formation of black holes with masses less than a solar mass. 

A black hole will unavoidably "devour" any matter in the immediate vicinity of its gravitational 
radius (more precisely, within about 3r g ). Therefore, if a PBH is located inside a star, the life of the 
star will be shortened by some time. After the passage of this time, all the remaining matter of the 
star will unavoidably collapse into the black hole. 

One consequence of the collapse of the star into such a PBH should be a powerful gamma-ray 
outburst, which is emitted by the remnants of the stellar material. It is possible that such gamma-ray 
bursts contribute to all the cosmological gamma-ray bursts that are observed in the Universe on a 
continual basis. 

Knowing the rate of accretion of matter onto a PBH and the distribution of PBHs in space, we 
can estimate the rate of associated gamma-ray bursts and compare it with the observed rate, and also 
estimate the number of PBHs in a galaxy. 



2 RATE OF ACCRETION OF A STAR ONTO A PBH 

A model for accretion onto a black hole has been considered in detail in pQ, [2], and we present here 
only the main conclusions of this theory. 

We consider the hydro dynamical theory of accretion onto a black hole. The gravitational radius r g 
of the PBH with mass m is assumed to be much smaller than the scale on which the self-gravitation of 
the star becomes stronger than the gravitation of the PBH. We denote this scale R sg ; note that it should 
be much smaller than the size of the star itself R@ ~ 7 • 10 10 sm, (for a star of mass M & ~ 2 • 10 33 g). 
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We then have: 

r g <^R sg <^R G , R sg = R Q lj?-J , r g = — — , R g = — ^ . (1) 

The rate of accretion onto a PBH with mass m is equal to (the Bondi solution): 

_ 47rA(Gm) 2 p Q 
(a 2 + PI 



c 3 / ,,2 j_ ;2 

Here, A is a constant of order unity, ca is the sound speed in the star (c 2 a 2 = dP/dp®), c(3 p is the speed 
with which the PBH moves through the star (at the perigee of its orbit, a « (5 P ), p G is the mean 
density of matter in the star. Thus, during its passage through the star, the PBH "digs a tunnel" in 
the stellar material with a cross-sectional area 

\irr 2 g ^ tttJ 

0p{a 2 + £2)3/2 ~ £4 ■ W 

The associated change in the star's mass is 

Am « p o sii » m • ^ g . (4) 

During this accretion, the internal parameters of the star (dependence of the temperature and density 
on radius) can be taken during computations to be equal to the solar values from [3]. Questions related 
to the stationarity of the accretion are considered in Appendix 1. 

3 CAPTURE OF A PBH IN A STELLAR ORBIT (ORBITAL 
CAPTURE) 

Moving in a galaxy with virial speeds of the order of c/3oo ~ 300 km/s (/?oo ~ 10 -3 ), a PBH has some 
probability to collide with stars. The cross section for such a collision S is it multiplied by the square 
of the maximum impact parameter at infinity, for which the PBH passes "within half the radius of 
the star. The quantity S is determined by equating half the stellar radius and the distance from the 
star to the PBH at the perigee of its orbit for the motion of the PBH along a hyperbolic trajectory: 

S = nR 2 Q ( \ + — %A « 2AkR% . (5) 



®U^2R Q (3l)~—"° 

Before the collision, the speed of the PBH is more than a third of the cosmic velocity for such a star. 
The speed at infinity (3oc is related to the speed just before the collision (at perigee) by the formula 

(3 2 p = Pl + § L , - &~ 2.3&0. (6) 

Since momentum is conserved in collisions involving the PBH, its speed will decrease as its mass 
increases. Thus, its speed immediately after a collision will be lower than it was just before the collision. 
The collisions will continue until the PBH is captured into an orbit by a star, i.e., until the speed of the 



PBH turns out to be less than a third of the cosmic speed after a collision: = JR 9 /Rq ~ 2.1^^, 
Thus, the total change in the speed of the PBH required for orbital capture should be 

A/3 = P - h « 0.2/Joo . (7) 
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In accordance with the conservation of momentum and expression the change in speed during the 
collision process is determined by the formula 



Am 3r„R n r a R 



A&«a, m "« « 16 ^ . (8) 



Accordingly, the number of collisions required for orbital capture is 

We can see that, in order for an orbital capture to occur as a result of the first collision, the mass of 
the PBH before the collisionmust m ~ 0.2M Q . In this case, the mass m of the PBH will approximately 
double result of the collision. 

4 RATE OF ORBITAL CAPTURES OF PBH IN A GALAXY 

We can determine the mean time T\ between collisions using the formula S'/3 00 cTi = Vb, where Vq = Rq 
is the mean volume per star in the galaxy (Rq ~ 10 light years ~ 10 19 sm): 

T x = « 3.8 • 10 23 years . (10) 

Thus, the time for the mean free path of a PBH (until its orbital capture) is r = k -T\. If the total 
number of uncaptured PBHs in the galaxy is N^h, some number of them (Nt) will be captured by 
stars over a time t: 

N t = N bh • [1 - exp(-t/r)] « N bh ■ t/r . (11) 

The number of PBHs in the galaxy is given by the relation Nbh = nMtat/fn, where Mtot ~ 10 46 g is the 
total mass of the galaxy (assuming it is similar to the MilkyWay) and k is a coefficient determining 
the mass fraction of PBHs in the galaxy (0 < k < 1). We can see that (jlip does not depend on the 
mass of the PBH. The rate of orbital captures of PBHs in the galaxy is then 

dN t _ s T 2AkcR 9 2GM tot 1 

dt ~ * ~ VqPI ' c 2 ~ K ' 10 6 years ' [ ' 

Thus, the rate of orbital capture can reach values of the order of one capture every million years in a 
MilkyWay-type galaxy. 

5 LIFETIME OF A PBH IN ORBIT AROUND A STAR 

If the lifetime of a PBH in orbit is less than the lifetime of the star (~ 10 9 yrs), PBHs will not be 
accumulated in the galaxy. 

The lifetime of a PBH in orbit is determined by the increase of its mass by a factor of two due 
to accretion. The corresponding number of revolutions n of the PBH around the star is given by the 
relation 



m _ 37PtR 2 & 2M 
n = — w — ^ ps . (13) 

Am Rg r g m 



The period of revolution about the star To is determined from the assumption that the major axis of 
the orbit is ~ 0.5i?o (half the mean distance between stars): 

?3/2 



r « irR%y(cJRg) « 1.7 • 10 8 yrs . (14) 
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We thus obtain for the lifetime in the orbit 

117/2* / R„ \ ( R \ 3/2 M o 

T = nr « g^J3 • M ■ ^ « 3.5 • • 10 8 years . 15 

\r g J \RgJ m 

The masses of PBHs can be constrained from above by data available from gravitational microlensing 
observations: m < O.IMq (see [6]). Therefore, in order for the time To to be less than the lifetime 
of the Universe (~ 15 • 10 9 yrs), the mass of a PBH must lie in the very narrow range of 0.03M Q to 
O.IMq. In this case, gamma-ray bursts induced by the collapse of a star into a PBH will occur from 
time to time in a galaxy. 

However, even in this case, this mechanism will not make an appreciable contribution to the 
observed gamma-ray bursts in a galaxy, since the rate of such events (j!2p is much less than the mean 
rate of gamma-ray bursts in a galaxy ~ 0.001 -j- 0.01 yrs~ x . 

It follows from the above analysis that PBHs can comprise any mass fraction of the dark matter 
in a galaxy. 



6 Appendixes 

6.1 TIME FOR ESTABLISHING A QUASI-STATIONARY ACCRETION REGIME 
AS A PBH IMPACTS A STAR 

The relative speed between the PBH and the particles in the stellar medium changes as a PBH 
moves through a star uniformly and linearly. Consequently, there are relative accelerations, which 
correspond to inertial forces. If we denote the coordinate of the PBH along its trajectory y = cf3 p t 
and the minimum distance between the particles in the stellar medium and the PBH trajectory x, this 
acceleration g can be expressed in terms of x and y as follows: 

I | _ jc(3 p x) . . 

191 ~ (X 2 + y 2 )3/ 2 • 1 j 

Here, it is convenient to assign the initial value of y to be at a point that is closest to a particle of the 
medium. The condition for capture of a particle in the medium by the passing PBH is given by the 
inequality 

J<TT2' (17) 
x A + y £ 

This corresponds to the dominance of the gravitational force over the inertial force, and must be 
satisfied at any moment in time. We thus obtain 

x < x ° = Tk ■ ( 18 ) 

This qualitative expression approximately corresponds to the cross section s = ttxq [see (|3|)]. 

The relaxation time for the stellar material during a collision with a PBH is determined (to order 
of magnitude) by the free-fall time for a particle to fall from a distance xq onto the PBH: 

"-safer (19) 

We have for a mass of m < O.IMq the time to <~ 10s. Thus, this time is much shorter than the time for 
the passage of the PBH through the internal layers of the star: [~ R@/(cf3 p )] ~ 1000s. Consequently, 
the accretion process very rapidly becomes quasi-stationary. 
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6.2 STATISTICS OF THE MASS AND SPEED DISTRIBUTIONS OF PBHS IN 
THE UNIVERSE 

The nature of PBHs in the early Universe (if indeed they exist) is not clear, and there are also no precise 
models describing their statistical distribution. Moreover, we cannot consider the early Universe to 
be in some kind of thermodynamic equilibrium. Therefore, the arguments presented below do not 
pretend to be the only valid description of the statistics of PBHs, and this has led us to present them 
in a separate appendix. 

Let us consider a system of a larger number of PBHs in a spatial volume V. We will consider this 
volume to be sufficiently large so that we can apply a statistical analysis to the PBHs. Their masses 
are bounded from below by the relation m m i n = E m i n /c 2 because of the condition for the Hawking 
evaporation of black holes. 

Let us assume that the PBHs have a canonical Gibbs distribution [5]. Thus, the probability density 
for the existence of a PBH with a mass m and a spatial velocity v is 

- = c 1 exp|- 7 ^=}, (20) 

where y = mc 2 / E m i n is the dimensionless rest energy of a PBH, a = E m i n /T is the dimensionless 
inverse temperature of the system, (3 = v/c the speed of the PBHs in units of the speed of light, and 
ci is a normalization constant. 
The normalization condition is 

oo 1 V 

J dy J 4tt(3 2 d(3w J d 3 r = 1. (21) 
10 o 

We now introduce a number of functions that depend on a: 

1 , ^ 



I n {a) = J 4irp 2 dp y/l-P 



o 

The condition (1211) then takes the form 



eXp (~ ^i% 2 ) ' I» = ~In-i{oc). (22) 



c x Vh = a. (23) 
We now introduce the mean energy of the PBH: 

U=( , V J)=7dy [ 47rff 2 d/3 -=y=c 1 Vexp { ■ (24) 

After integrating over y, this expression can be written in the form 

Ua 2 = Cl V [ha + Ji] . (25) 

The energy and entropy of any system are additive quantities, equal to the sums of the contributions 
from their individual parts. Black holes should not be an exception to this behavior. Therefore, the 
total energy and entropy of the system will have the forms 

U tot = Jdy [ 4vr/? 2 d/3iV^=L==c 1 Uexp|- ^__ | = NU . (26) 



o 

j 

Stot 



( ~1 

Jdy J 47T/3 2 d/?iYS oCl y exp <j - y==f f = N ^ ■ ^ 
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where N is the total number of PBHs in the system, So = ciy 1 is the entropy of a PBH with mass m, 
and C2 is a dimensionless coefficient [3J: 

C2 = » i at mmin » i -5 gm 28 

nc 

Thus, 5tof = (So) Utot/U . Using the thermodynamic expression for the temperature, 1/T = dStot/dE to t 
(see [5]), we obtain 

oo l , >. 

a = — = jdyj4nP dP—y expj-^—j. (29) 
10 *■ 7 

Integrating over y, we obtain 

[/a 4 = C2Cl V [2/3 + 2a/ 2 + n 2 h} . (30) 



Equations ([23]) . ([25]) and ([30]) can be used to find a and U, e.g., by multiplying ([25]) by a 2 and 
substituting this into (|30p . We differentiate the resulting expression with respect to a to obtain 

a 2 /_x - a(c 2 + 2)I - 2h = (31) 

This is a transcendental equation in a. However, if C2 3> 1, then a> 1. We will use this fact and 
multiply the function I n by the factor exp{a} /(47r). Denote the new functions I n : 



I n (a) = I p 2 d/3 



■mi/ 1 \ 1 / 02 



1-/3 



exp ja II - ^=f j I « exp{-a/3 2 /2} • /3 3 • h - ^ j . (32) 



In 7 n , the integral is already rapidly cut off by the exponential when /3 2 ~ 2k/ a <C 1, where k is a 
coefficient of order 10. Thus, the functions obtained can be expanded in a series in /3q, keeping only 
the first term. The function I n preserves the form of (|3ip . Substituting the asymptotic (|32p into (|3ip 
yields 

a 2 -a(c 2 + 2-K)-2 + 2K/a = 0. (33) 
We find in the main approximation 

a « c 2 , 17 » 1 + 5/(4a) » 1 . (34) 

Hence, the temperature of the Gibbs distribution of PBHs is of order the Hawking temperature (T H ) 
for m min : 

T =l^^Ejran = ^ = = ^ ^ / lO 2 ^ \ Q 

a c 2 47rGm m j n \m min J 

It follows that the bulk of PBHs have very narrow ranges of their mass (from m m i n to m m i n {l + 1/a)) 
and velocity: 

10 — ^ 

from to A) < /3oo < 1. 

The probability density for the PBHs falls off exponentially outside these intervals. With such a low 
value for (3q, its contribution to the motion of the PBH can be neglected compared to the virial speeds 
in the galaxy, which are determined by the gravitation of all its bodies. 

Thus, to a high degree of accuracy, this model for PBHs represents an ideal gas with particles of 
mass m m i n and velocities of the order of the virial velocity in the galaxy. 
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